Relaxation timescales and decay of correlations in a long-range interacting quantum 

simulator 
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We study the time evolution of correlation functions in long-range interacting quantum Ising 
models. For a class of initial conditions, exact analytic results are obtained in arbitrary lattice 
dimension, both for ferromagnetic and antiferromagnetic coupling, and hence also in the presence 
of geometric frustration. In contrast to the nearest-neighbor case, we find that correlations decay in 
time like stretched or compressed exponentials. Provided the long-range character of the interactions 
is sufficiently strong, pronounced prethermalization plateaus are observed and relaxation time scales 
are widely separated. Specializing to a triangular lattice in two spatial dimensions, we propose to 
utilize these results for benchmarking of a recently developed ion-trap based quantum simulator. 
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When a physical system is coupled to a heat bath, one 
expects to observe thcrmalization to an equilibrium state 
whose temperature is determined by the bath properties. 
For an isolated many-body system, i.e., in the absence 
of a heat bath, the situation is less clear, although some 
kind of relaxation to equilibrium may be expected for suf- 
ficiently large generic systems and suitable observables. 
Recent progress in experiments with cold atoms and ions 
[1] has stimulated intense theoretical interest in equilibra- 
tion and thermalization behavior of isolated many-body 
quantum systems [2, 3]. General mechanisms leading to 
thermalization have been proposed [4] , rigorous proofs of 
equilibration have been obtained for generic Hamiltoni- 
ans [5], and analytic as well as numeric model studies 
have been reported (see [3] for a list of references). Much 
less is known about the time scales on which relaxation 
to equilibrium takes place. In this Letter we study the 
time evolution of correlation functions in isolated long- 
range interacting quantum Ising systems on arbitrary lat- 
tices. We derive exact analytic results that further our 
understanding of different relaxation timescales and gen- 
eral noncquilibrium properties of long-range interacting 
systems. 

Apart from broadening the theoretical understanding, 
our results may also prove beneficial for the interpreta- 
tion of data from cold atom or ion experiments. Special- 
izing to a triangular lattice in two spatial dimensions, we 
propose to utilize our results for benchmarking a recently 
developed ion-trap based quantum simulator [6]. This 
quantum simulator consists of several hundred beryllium 
ions stored in a Penning trap and confined to a planar 
slab. Due to their mutual electrostatic repulsion, the ions 
arrange into a two-dimensional Coulomb crystal on a tri- 
angular lattice (see Fig. 1). The valence electron spins of 
the 9 Be + ions are the relevant degrees of freedom used 
for quantum simulation. Effective interactions between 
these spins are induced by transverse motional modes of 



the Coulomb crystal. Under the assumption of small, 
coherent ion displacements, it was shown in Refs. [7, 8] 
that the resulting interactions are described by the Ising 
Hamiltonian 



H 



i<j 



B„ 



(i) 



Here i and j label the N ions on the triangular lattice, 
(Ti = (erf, erf, of) denotes the vector of Pauli matrices for 
ion i, the Jij arc coupling coefficients, and is an ef- 
fective magnetic field generated by externally applied mi- 
crowave radiation. Unlike in the conventional Ising model 
[9], spin-spin coupling is not restricted to nearest neigh- 
bours on the lattice, but extends over all pairs of ions. 
The coupling coefficients Jij can be expressed in terms 
of the transverse phonon cigenfunctions of the lattice and 
a few other experimental parameters (Eq. (4) of [10]). A 
numerical evaluation of that expression for the given lat- 
tice geometry shows that the approximation Jij oc D~" 
holds to a very good degree [6] . Here Dij denotes the Eu- 
clidean distance between sites i and j on the lattice, and 
a is an exponent that can in principle be tuned within 
the range < a < 3 by varying the difference in fre- 
quency of two off- resonant lasers used in the experiment. 
The absolute values and even the signs of the coefficients 
Jij can also be tuned, allowing for the investigation of 
ferromagnetic as well as anti-ferromagnetic couplings. In 
the latter case, due to geometric frustration on the tri- 
angular lattice, spin liquids and other exotic quantum 
phases may possibly occur. It should be feasible to ef- 
ficiently simulate the behavior of a long-range quantum 
Ising model on a triangular lattice with several hundred 
trapped ions, well beyond system sizes that can be dealt 
with in numerical simulations. However, this will be true 
only if the assumptions used in [8] for the derivation of 
(1) are indeed satisfied in the experiment. To date, the 
Penning ion trap simulator has been benchmarked only 



through comparison with mean field theory, a classical 
limit of (1) where quantum fluctuations and correlations 
are ignored. The exact results of this Letter enable a 
much higher level of benchmarking, and a determination 
of the time scales over which the expected emulation of 
quantum effects of an Ising model is indeed realized. The 
results can also be applied to the linear radio-frequency 
ion trap quantum simulators [11] and, along similar lines, 
benchmarking of quantum magnets emulated by means 
of ultracold molecules should be possible [12]. 

For the theoretical analysis of relaxation times and the 
decay of correlations, we study a model which is more 
general than the Hamiltonian in (1) in some respects, 
and more restricted in others: We generalize to lattices of 
arbitrary spatial dimension d and arbitrary lattice struc- 
ture, but come back to the two-dimensional triangular 
lattice later in this Letter. The effective magnetic field 
B^ in (1) is required to point in z direction, yielding a 
Hamiltonian of the form 

i<j i 

where the index £ indicates the presence of a longitudi- 
nal field B. At first sight it may look as if the resulting 
problem is purely classical, as all terms in H? commute 
with each other, and equilibrium properties such as the 
partition function are indeed identical to those of the cor- 
responding classical Ising model. For noncquilibrium cal- 
culations, however, quantumncss enters through obscrv- 
ables and/or initial density operators that do not com- 
mute with the Hamiltonian (2), and entanglement and 
other genuine quantum properties may be present. Ob- 
scrvablcs consisting only of of operators commute with 
the Hamiltonian and therefore show no relaxation behav- 
ior, but other obscrvablcs do, as we will confirm in the 
following. In spite of its simplicity, the model (2) exhibits 
non-trivial and remarkably rich nonequilibrium behavior. 

To study the relaxation to equilibrium in this model, 
we want to compute the time evolution of spin-spin cor- 
relations 

(ofof) (t) - Tr (e^of of e- iff «%) (3) 

in units of h = 1. A significant simplification is achieved 
by restricting the initial states to density operators po 
that are diagonal matrices in the a x tensor-product 
eigenbasis. This idea has been used in [13, 14] for the 
computation of expectation values of a rather restricted 
class of 1-spin observables. In the present Letter, these 
restrictions are relaxed and the exact analytic calcula- 
tions are extended to correlation functions for arbitrary 
finite system sizes. Details of the calculation are reported 
in the Supplementary Material. For vanishing longitudi- 
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FIG. 1. (Color online) Left: A top-view resonance fluores- 
cence image showing the center region of an ion crystal cap- 
tured in the ions' rest frame. Fluorescence is an indication of a 
valence electron being in the spin-up state. Here, all ions are 
in the spin-up state. The lattice constant is approximately 
20 jj,m. Center and right: Triangular lattices on hexagonal 
patches of side lengths L — 8 and 16. 

nal magnetic field B = 0, the final results are 

<«) W = <ofo?>(*) = <*?o?)(t) = 0, (4a) 
(ofaf)(i) = (of)(0)sin(2*J i>i ) JJ cos(2tJ ilfc ) , (4b) 

(afa*)(t)=Pr j+ P+, (^)(t) = i% - P+, (4c) 
with 

p t, = \ °?>(o) n c ° s t 2 ^ ± *i • ( 4d ) 

These expressions are valid for arbitrary coupling con- 
stants Jij, and therefore apply in arbitrary spatial di- 
mension and on arbitrary lattices. Moreover they can 
easily be evaluated on a personal computer for millions 
of spins. The presence of a longitudinal magnetic field 
B^O modifies Eqs. (4a)-(4d) by imposing an additional 
spin precession at an angular frequency proportional to 
B (Supplementary Material). For the relaxation to equi- 
librium we are interested in, precession turns out to be 
irrelevant and we will discuss the B = case in the fol- 
lowing. A generalization to arbitrary n-spin correlation 
functions is also straightforward. Since any operator on 
(D 2 can be expanded in terms of products of Pauli op- 
erators, such results for correlation functions of arbitrary 
order permit in principle the computation of the time- 
evolution of arbitrary operators. This docs not add much 
to the discussion of the physical phenomena we are inter- 
ested in here, and we will therefore restrict the discussion 
to 2-spin correlations. 

As an example, the time evolution of spin-spin cor- 
relations is illustrated for system parameters similar to 
those of the ion trap experiments of Britton et al. [6] : We 
consider hexagonal patches of triangular lattices (Fig. 1), 
and Jij = JD^f proportional to the crth power of the 
inverse Euclidean distance of sites i and j on that lat- 
tice, with J being a coupling constant. The plots in Fig. 
2 suggest that correlations decay to their microcanonical 
equilibrium values (afa^) mc = where a,b £ {x, y, z}, 
but this decay is only apparent: Since the expressions in 
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FIG. 2. (Color online) Time evolution of the (normalized) 
spin-spin correlations as given by (4a)-(4d). The lattice sites 
i and j are chosen one lattice site to the right, respectively left, 
of the center of the hexagonal patch, as indicated by the blue 
dots in Fig. 1 (center and right). The various curves of the 
same color in each plot correspond to side lengths L = 4, 8, 16, 
32, and 64 (from right to left) of the hexagonal patches. The 
left panel is for power law interactions with exponent a = 1/2, 
but results are qualitatively similar for all a between zero and 
d/2. The right panel is for a = 3/2, with qualitatively similar 
results for all a > d/2. The unit of time is 1/J. 



(4b)-(4d) consist of products of TV trigonometric func- 
tions, the evolution in time of spin-spin correlations is 
quasi-periodic for any finite number TV of spins, and re- 
laxation to equilibrium cannot occur in the strong sense. 
Instead, recurrences will repeatedly bring the system ar- 
bitrarily close to its initial state. These recurrences occur 
on a time scale that is exponentially large in TV and, al- 
ready for moderate system sizes, do not show up on the 
time scales plotted in Fig. 2. (For example, for a lattice 
with side length L = 4 and a = 3/2, the first significant 
recurrence occurs at a time three orders of magnitude 
larger than the relaxation time). For a mathematical 
analysis, however, the recurrences may pose an obstacle. 
In order to study relaxation to equilibrium in a meaning- 
ful way, one can either resort to a weaker (probabilistic) 
notion of equilibrium [5], or, and this is the strategy we 
use, suppress the quasi-periodic recurrences by working 
in the thermodynamic limit TV —> oo. With (4b)-(4d) as 
a starting point, we have derived upper bounds on the 
time evolution of spin-spin correlations in the thermody- 
namic limit for a few different lattice structures, and it 
turns out that they are all of the form 



<afaj>(t), <afajf>(i)<Pr 



with 



P; 



ho 



c )(0)exp[-C± (a)N q±( > a h p±( -^ 



(5a) 



(5b) 



with q + (a) = max{0, 1 — 2a/d}, q~(a) = 0, p + (a) = 
min{2, d/a}, and p~(a) = d/(l + a). The constants 
C depend on the exponent a and the lattice structure. 
Other 2-spin correlation functions have similar proper- 
ties, and details of the derivation are provided in the 
Supplementary Material. For the example of a trian- 
gular lattice in two dimensions with a = 3/2 and i,j 
as in Fig. 1, we obtain C+- « 11.04 and Cr. k, 1.119, 
and the corresponding bound (5a) is compared to exact 



finite-system results in Fig. 3 (left). In contrast with the 
exponential decay of correlations in the nearest-neighbor 
Ising model [15], the bound (5b) reveals a compressed 
exponential decay for a < d/2, and a stretched exponen- 
tial decay otherwise — even for large a, where power-law 
couplings are expected to correspond to nearest-neighbor 
interactions in the limit a — > oo. 

There are other, and even more striking, a dependent 
aspects of the relaxation dynamics: The plot in the right 
panel of Fig. 2, representative for exponents a > 1, shows 
a simple relaxation to equilibrium with a single relevant 
time scale. The left panel of Fig. 2, representative for ex- 
ponents < a < 1, differs considerably: In a first step, 
the <rf <rj and cfcrj correlations display a Gaussian decay 
on a fast timescale n to precisely 1/2 of the afaj initial 
value, before finally relaxing to their vanishing equilib- 
rium value on a much longer timescale r%. t\ coincides 
with the timescale on which 1-spin observables (green 
lines in Fig. 2) relax, while correlations are not yet equi- 
librated. Then the system remains prethermalized [16] 
for a relatively long period of time, a behavior observed 
in the left panel of Fig. 2 as a conspicuous plateau. Dif- 
ferent from earlier observations of prethermalization in 
quantum dynamics, the ratio t\/t2 goes to zero in the 
large- TV limit, i.e., separation of timescales become more 
pronounced with increasing system size. We believe that 
such a behavior is peculiar to long-range interacting sys- 
tems, as it bears some similarity to divergent (with TV) 
relaxation timescales as observed for classical long-range 
systems in the astrophysical context [17]. Similarly, the 
separation of timescales becomes more pronounced also 
in the limit a — > 0, in which the slow relaxation timescale 
diverges. This fact may be explained by the presence of 
TV(TV— 1) additional symmetries in the a — case, where 
all operators of oj+of crj commute with the Hamiltonian 
(2). Although these symmetries are absent for a > 0, the 
slow timescale can be seen as a remnant of the weakly de- 
stroyed a = symmetry. While 2-spin correlations may 
exhibit relaxation in a two-step process, we have checked 
that no third timescale emerges when the calculations are 
generalized to 3-spin correlations. 

Multiple timescales of relaxation may emerge for var- 
ious reasons. One possible scenario that has been ob- 
served previously, both in theory and experiment [18], is 
that dephasing is responsible for prethermalization, while 
a collisional mechanism causes relaxation on a slower 
timescale. In order to test which of these mechanism is 
at work in the long-range Ising model, we have computed 
the time evolution of the n-spin purity 



7n(*) = Tr (pi (t)) 



(6) 



where p n is the n-spin reduced density operator, as ob- 
tained by tracing the density operator p over all but n 
spins. Knowledge of the 1-spin expectation values (erf), 
(of), and (erf) allows for the reconstruction of p\, and 
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FIG. 3. (Color online) Left: Upper bound (5a) (straight red 
line) compared to the exact time evolution (4c) (wiggly blue 
line) of the spin-spin correlator (afaj}(t)/(afaj}(0), plotted 
vs. rescaled time t 2 ^ 1+a K The linearly decaying trend in the 
plot, superimposed by fluctuations, confirms the asymptotic 
behavior over a range of more than hundred orders of mag- 
nitude. Right: Time evolution of the 1- and 2-spin purities 
71 and 72 as defined in (6). Each relaxation step in spin-spin 
correlations in Fig. 2 is accompanied by a drop in 72. This 
suggests that dephasing is responsible for all relaxation steps. 

additional knowledge of the spin-spin correlations (4a)- 
(4d) facilitates reconstruction of p2- For the long-range 
Ising model, as shown in Fig. 3 (right) for a = 1/2, both 
relaxation steps of spin-spin correlations turn out to be 
associated with a drop in the purity 72 , and this is an in- 
dication that both relaxation steps are caused by dephas- 
ing. In contrast, a collisional relaxation would be charac- 
terized by a redistribution of mode occupation numbers 
or probabilities and requires inelastic scattering processes 
between these modes. Since the Ising model (2) has only 
commuting terms, no inelastic collisional mechanism is 
available. 

In the trapped ion simulator of [6], spontaneous emis- 
sion from the off-resonant laser beams used to engineer 
the long-range Ising interaction (1) sets one experimental 
time scale that, over longer intervals of time, the system 
can no longer be considered isolated. However, we cal- 
culate that with upgrades to the set-up described in [6], 
the quantum many-body relaxation time scales can be 
short compared to the spontaneous emission time scale. 
Specifically, with N = 217 ions and an optical dipolc 
force generated from two 10 mW beams crossing with an 
angular separation of 35° (see Fig. 1 of [6] ) and frequency 
difference tuned to generate an a = 1/2 power law inter- 
action, we calculate the two relaxation timescales cor- 
responding to P^j in (4d) to be ~ 30 fjs and « 430 /xs, 
respectively. The spontaneous emission time for this con- 
figuration is (1/r) ~ 4 ms. All of these time scales are 
short compared to the T2 > 50 ms coherence time of the 
Be + valence electron spin qubits. 

To compare with the exact results reported here re- 
quires an experimental measurement of the spin-spin cor- 
relations after the Ising interaction has been applied for 
a variable period. This is readily accomplished with 
trapped ions. Reference [6] uses state-dependent res- 
onance fluorescence to measure spin orientation in the 
<j z basis. Specifically global fluorescence collected from 



all the ions in the array is used to measure the global 
spin state of the system (i.e., the total number of spins 
in the |f) z state and the total number of spins in the 
14,) state). Shot-to-shot fluctuations in these measure- 
ments are sensitive to the second-order moment (J z ) of 

the total ^-component of the spin J z = J2iLi °f/2- Of 
particular interest are measurements of the second-order 
moments in directions perpendicular to the mean com- 
posite spin vector. For example, with all spins initially 
pointing along the x-axis, fluctuations in measurements 
of J z after rotation about the x-axis by an angle 9 are 
sensitive to 

N 1 

- sin 9 cos 9 (a^ a] + <t?<t?) + cos 2 9 (ofo*)] . (7) 

These spin-squeezing measurements [19] are sensitive to 
pairwise correlations summed over all the spins in the en- 
semble. In addition to global fluorescence measurements, 
the top-view image in Fig. 1 (left) indicates that individ- 
ual spin detection may also be possible in this system. 
Individual spin detection is complicated by the rotational 
stability of the crystal, but has significant advantages in 
the reduction of potential sources of noise in the second- 
order moment measurements, as well as permitting a mi- 
croscopic measurement of the spin state of the system. 

Once such a quantum simulator is benchmarked, it 
may in turn be used for the testing of more general 
approximate theories. Moreover, there are certain phe- 
nomena which are believed to be peculiar to long-range 
interacting systems and which might receive their first 
experimental confirmation in the ion trap set-up. One 
example is the threshold at a = d/2 below which a sec- 
ond timescale of relaxation emerges. There is evidence 
that not only the existence of this dynamical long-range 
threshold, but also its numerical value d/2, is universal 
for classical as well as quantum mechanical long-range 
interacting systems [20]. Other long-range peculiarities 
include nonequivalent equilibrium statistical ensembles, 
a phenomenon that has been known in the astrophysical 
context for decades [21] but has not seen experimental 
verification in a tabletop experiment. 

In summary, under convenient restrictions on the ini- 
tial conditions, we have obtained exact analytic results 
for spin-spin correlations of long-range quantum Ising 
models in arbitrary spatial dimensions and for various 
lattice structures. Prethermalization, widely separated 
timescales, and other phenomena that further our un- 
derstanding of the timescales governing the relaxation 
to equilibrium were discussed. Comparison of our ana- 
lytic results to experimental measurements of spin-spin 
correlations will enable benchmarking of the trapped-ion 
quantum simulator of [6] in a regime where quantum 
effects are paramount. Subsequently, the benchmarked 
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quantum simulator may be employed to investigate re- 
laxation timcscalcs in parameter regimes and for initial 
states where analytic results presently are not available. 
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Supplementary Material 



Appendix A: Exact calculation of spin— spin 
correlations for finite N 

Starting point of the calculation is the general expression 
of the expectation value 

(of erf >(*) = Tr (e^Vf of e^Vo) (A.l) 

with respect to the initial density operator po, where 
we assume po to be diagonal in the a x tensor prod- 
uct eigenbasis. Spin ladder operators arc defined as 
a = (a x ±ia y )/2. Since all terms in the Hamiltonian He 
commute, we can factorize the time evolution operator, 

exp (-iH e t) = Y[ exp (iJ fc ,;o|of i) JJ exp {\Ba z m t) , 

k<l m 

(A.2) 

and similarly for the Hcrmitian conjugate. All factors 
in (A.2) that contain neither af nor erj commute with 
a i o" . . To compute the time evolution in (A.l), we there- 
fore have to deal with the expression 

[ exp (-iJ^o-f ofi) exp (-Ukja^ap) exp (-\Ba z t) 

x exp (—\Ba z t) erf erf exp (iBojt) exp (iBaft) 

x JJ exp (iJijcrf aft) exp (i J M of aft) . (A.3) 

Making use of [a z ,a ± ] = ±2cr ± , the time evolution due 
to the magnetic field B simplifies to 

exp (— LBof t) exp (—\Baft) afaf exp (iBa z t) 

x exp (i.Bof t) = erf erf exp (T&Bt) . (A.4) 

Picking one lattice site k ^ the time evolution of 
of of due to the interaction with the spin at k can be 
written as 

exp (-iJ kii a%a?t) exp (~U k ^a z k a z t) afa j 
x exp (\J k>j a%ap) exp {\J k>i a%aft) 
= afaf cos [2t(J hk + Jj. k )} 

T iofofok sin [2t(.7»,* + J J,fc)l • ( A - 5 ) 

Since the initial state po is assumed to be diagonal in 
the a x tensor product eigenbasis, only diagonal elements 
of the operator e lHet af of e~ lHet in the same basis con- 
tribute to the trace in (A.l). For this reason we can drop 
the second term on the right-hand side of (A. 5), as it is 
proportional to af.. Inserting (A.4) and (A. 5) into (A.l), 
we obtain 

(°t<rf >(*) = (ofof >(0)exp( T 4LBt) 

x ] [ cos[2t(J i:k + J j>k )] , (A.6a) 



where (of of )(0) = Tr (of of po). A similar calculation 
yields 

(°t°7)® = n c ° s tw.* - j ^ > 

Mm 

(A.6b) 

a notable difference to (A. 6a) being that the presence of 
an external magnetic field B has no effect on this expec- 
tation value. From 

<^:i«::o>;. (A.7a) 
4cf of = of a* + afaj T iof (7? ± iafaj, (A.7b) 

we obtain 

afa] = 2$t(afaJ+afaf) ) (A.8a) 

m{o?vf-otof), (A.8b) 

of of = T23 (of of - trf o-f ) , (A.8c) 

erf e^ = ±23 (afaf + afaf) , (A.8d) 

and furthermore, using the properties of the initial den- 
sity operator po, 

i(afaf)(0)=A(afaf)(0) = (afaf)(0). (A.9) 

Inserting (A.6a), (A.6b) and (A.9) into (A.8a)-(A.8d), 
we arrive at the final expressions 

(afaf) (t) = P~ + COS (ABt)P+ , (A.lOa) 
<erferj)(t) = Pr. co S (4Bi)P+, (A.lOb) 
(afa y j )(t) = -sm(4Bt)P+, (AlOc) 

with Pf as defined in (5d). A similar calculate yields 

(a*a z )(t) = S m(2Bt)Pi,j, ("MX*) = ^s(2Bt)P z v 

(A.lOd) 

with 

Pi,j = ~ (of >(0) sin(2y y ) [] cos(2tJ iifc ). (A.lOe) 

Finally, from symmetry considerations we obtain 

(<r?erf )(i) = 0. (AlOf) 

A generalization of these calculations to n-spin corre- 
lation functions is straightforward, yielding for example 

(afafaf)(t) = (af af of >(0) exp ( T 6iP4) 

x [J cos[2t(J itl + J jt i + J kt i)], (A.ll) 

and similar expressions for other correlation functions. 
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Appendix B: Upper bounds on spin— spin 
correlations in the thermodynamic limit 

In the thermodynamic limit of infinite system size, a 
bound on the product 



JV 



i -J 



[ | cos [2 (J fc ,i± J k ,j)i\\ 



(B.12) 



k=l 



in Eq. (5d) with Jjj = D i Q and a > is derived. For 
any given t, one can find a compact region g i -(t) (con- 
taining i and j) of the infinite triangular lattice G such 
that 

|2(J M ± J kd )t\ < tt/2 (B.13) 
for all k G G \ gfj{t). Since | cosx| < 1, we have 

II |cos[2(J fc)i ±J fc)J -)t]|. (B.14) 



Next, using the inequality 

2 

< exp 



,2x 

cosx\ < 1 — I — 

7T 



for x < 7r/2, 
(B.15) 



Each individual summand (Jk,i ± Jk,j) can be minorized 
by (Jftr.i ± Jk.j) 2 , where K is determined from k by go- 
ing vertically upwards along the dashed (red) line in the 
above illustration until a boundary point of the shaded 
region is reached. To a boundary point labeled by K(r) 
there correspond r + 1 points in the sum, and this allows 
us to bound the sum by 



— a 
3,K(r 



E (Jk,i±Jk, j ) 2 >2J 2 J2 r ( D 7,K(r) ±D : 
keG\gf-(t) r=R±(t) 

(B.18) 

To exclude lattice sites in the region gfj{t), i.e., in order 
to satisfy the inequality (B.13) used earlier, R${t) has 
to be chosen large enough. An asymptotic analysis of 
(B.13) shows that this condition can be implemented by 
choosing 



R+{t) > 



1/q 



„ , , { AaJt\ 



(B.19) 

where the condition on R${t) works always, the one for 
Rq (t) only for sufficiently large t. Inserting the distances 



we obtain 



D itK{r) = r, Dj, K {r) = V5 2 +r5 + r* (B.20) 
and bounding the sum by an integral, we obtain 



Vf, < exp 



Wat, jit) 



(B.16) 



For simplicity, we consider lattice sites i and j symmet- 
rically arranged to the right and left of lattice site 0, as il- 
lustrated below for the example of distance 8 = Di j = 2, 
but other arrangements can be treated in a similar way. 

Km 




Disregarding all sites that lie outside the grey shaded area 
hij in the above illustration, we obtain a lower bound on 
the sum in (B.16), 

E (Jk.i± Jk,j)' 2 > E (Jk,i± Jk,j) 2 ■ 

fceG\ 9l ± (t) keh z _ J n(G\ g f j (*)) 

(B.17) 



Vf tj < exp 



32J 2 < 2 



. 1 1 

dr r I — ± — , 

r° V^ 2 + rS + r 2 



(B.21) 

valid in the limit of large R. Since we are interested in 
the limit of large R and t, we can expand the integrand 
to leading order in 1/r, obtaining 



1 



1 



V<5 2 + rS + r 2 
1 

\/5 2 +rS + r 2 ' 



4r 



(B.22a) 



(B.22b) 



Inserting these expressions into (B.21), the integral can 
be solved by elementary means. Making use of the 
asymptotic equality R 2 ~ A/3 and performing the limit 
A — y oo, we obtain the final results 



exp 



exp 



64J 2 t 2 
Tr 2 {l-a) 



A 

2/a 



and 



V l3 < exp 



8Jt 
n 



P_ / WA 2/(1+a) 
4a \ 7r J 



for a < 1, 

for a > 1, 

(B.23a) 

(B.23b) 
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FIG. 4. Logarithmic plot of the products (top) and V i j 

(bottom) as a function of rescaled time t p < - a \ evaluated for 
a hexagonal patch of a triangular lattice with side length 
L = 16 for a = 3/4 (left) and a = 5/2 (right). The ex- 
ponents p + (a) = min{2, 2/a} and p~(a) = 2/(1 + a) of the 
time rescaling are those predicted by the bounds (B.23a) and 
(B.23b) for the asymptotic behaviour at large L and t. The 
linear decaying trend in the plot, superimposed by fluctua- 
tions, confirms over a range of more than hundred orders of 
magnitude that the asymptotic behaviour of Vfj is indeed 
as given by the bounds. The straight lines in the plots are 
the bounds (B.23a) and (B.23b), indicating that the numeri- 
cal constants in the exponents of (B.23a) and (B.23b) are, as 
expected, underestimated. 



valid in the limit of large R and t. A comparison of 
these bounds with an exact evaluation of (B.12) for fi- 
nite lattices is shown, over more than hundred orders of 
magnitude, in Fig. 4. 



